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We investigate the electronic and magnetic properties of the frustrated triangular-lattice anti- 
ferromagnets CS2CUCI4 and Cs2CuBr4 in the framework of density functional theory. Analysis of 
the exchange couplings J and J' using the available X-ray structural data corroborates the values 
obtained from experimental results for Cs2CuBr4 but not for CS2CUCI4. In order to understand this 
discrepancy, we perform a detailed study of the effect of structural optimization on the exchange 
couplings of CS2CUCI4 employing different exchange-correlation functionals. We find that the ex- 
change couplings depend on rather subtle details of the structural optimization and that only when 
the insulating state (mediated through spin polarization) is present in the structural optimization, 
we do have good agreement between the calculated and the experimentally determined exchange 
couplings. Finally, we discuss the effect of interlayer couplings as well as longer-ranged couplings in 
both systems. 

PACS numbers: 71.15.Mb,71.20.-b,75.10.Dg,75.10.Jm 



I. INTRODUCTION 

For almost two decades the frustrated antiferromag- 
nets CS2CUCI4 and Cs2CuBr 4 have been considered as 
experimental realizations of a frustrated triangular lat- 
tice [T] [5] . Both systems crystallize in the space group 
Pnma [51 13] [see Fig. [I] (a)] and are characterized by a 
layered arrangement of Cu 2+ ions in a triangular pattern 
parallel to the be plane. The two-dimensional character 
of magnetic interactions between the spin-^ Cu 2+ ions 
was confirmed by neutron scattering and susceptibility 
measurements in both systems [SHS] and was successfully 
modeled [TU1 E] by a two-dimensional Heisenberg Hamil- 
tonian containing an anisotropic interaction term of the 
Dzyaloshinskii-Moriya type [12] [13] . 

In spite of their structural similarity, CS2CUCI4 
and Cs2CuBr 4 have rather different magnetic behavior. 
While in Cs2CuBr 4 magnetic excitations are localized 
and the field dependent magnetization exhibits two well- 
defined plateaux [5], CS2CUCI4 shows fractional spin exci- 
tations and spin liquid behavior over a broad temperature 
range, as revealed in inelastic neutron scattering exper- 
iments [5]. The dissimilar behavior between CS2CUCI4 
and Cs2CuBr4 has often been attributed to their unequal 
degree of frustration, determined as the ratio J'/J be- 
tween the inter-chain exchange coupling J' and the dom- 
inant intra-chain exchange coupling J in the underlying 
triangular lattice [see Figs, [I] (b) and (c)]. 

A ratio of J'/J = 0.74 [15] has been suggested for 
Cs2CuBr4 by comparing the ordering vector of a heli- 
cal incommensurate structure observed in neutron elas- 
tic scattering experiments with the one obtained from in- 
verse temperature series expansions for a spin-^ Heisen- 
berg model on an anisotropic triangular lattice [14] . In 
contrast, a ratio of J'/J = 0.34 [11] was derived for 
Cs 2 CuCl4 from comparison of spin- wave calculations for 



the spin-| Heisenberg model with the magnetic exci- 
tation spectrum observed in neutron scattering exper- 
iments in the presence of an external magnetic field 
far above saturation. These observations indicate that 
CS2CUCI4 is less frustrated and more one-dimensional 
than Cs2CuBr 4 . 

While a large amount of work has been devoted to the 
description of the fractional quantum states and to un- 
derstanding the phase transitions of the two-dimensional 
frustrated spin model with spatial anisotropy [Tf3H54] . a 
detailed comparative analysis of the electronic, magnetic 
and structural properties of these systems as well as a 
deep understanding of the origin of the different behav- 
ior is still missing. 

In this work, we present an extensive density func- 
tional theory (DFT) study of the microscopic proper- 
ties of CS2CUCI4 and Cs2CuBr4 and compare our results 
with experimental data. We consider different exchange- 
correlation functionals in order to also investigate the de- 
pendence of the electronic, magnetic and structural prop- 
erties on these choices. Our study of the performance 
of different exchange-correlation functionals is motivated 
by the fact that physical properties of recently discov- 
ered high temperature Fe-based superconductors are ex- 
tremely sensitive to the details of DFT calculations [3SJ- 
|4T] . Also, DFT studies of a recently topical layered Mott 
insulator TiOCl [42Tf5T] reveals that it is essential to use 
a suitable exchange-correlation functional in DFT calcu- 
lations to describe correctly the behavior of this system. 

Out of the electronic structure calculations we derive a 
tight-binding (TB) Hamiltonian and estimate the Heisen- 
berg exchange coupling constants between Cu ions from 
total energy calculations. Our DFT derived effective 
models incorporate a larger number of interacting Cu 
neighbors compared to the models used for the experi- 
mental data analysis. We show that some of these terms 
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are crucial for understanding the behavior of CS2CUCI4 
and Cs2CuBr 4 . 



II. STRUCTURAL AND ELECTRONIC 
PROPERTIES 

(a) Structural optimization schemes 

In this paper, we consider the experimentally measured 
structures of Cs 2 CuCi4 and Cs 2 CuBr 4 as well as internal 
coordinates obtained from the structural optimization 
with various exchange-correlation functionals. The ex- 
perimentally measured unit cell parameters of CS2CUCI4 
are taken from Ref.[3jand those of Cs2CuBr4 from Ref. [U 
Both structures were measured at room temperature. 

Crystal structure optimization of CS2CUCI4 and 
Cs 2 CuBr4 is required since calculations with the exper- 
imentally determined structural parameters find quite 
large interatomic forces, which indicates that the struc- 
tures are not in equilibrium within DFT. This holds 
for the calculations with all the considered exchange- 
correlation functionals listed below. In all cases, the 
interatomic forces are larger in Cs 2 CuCl 4 compared to 
Cs 2 CuBr 4 . 

Based on our structural optimizations, we will discuss 
in detail the electronic properties of both compounds. 
For the structure optimization we apply the following 
schemes with different approximations to the exchange- 
correlation functional within DFT and different magnetic 
configurations: (1) the spin-independent local density ap- 
proximation (LDA[nm]) [52] ; (2) the spin-independent 
generalized gradient approximation (GGA[nm]) [S3]; 
(3) the spin-dependent GGA with a ferromagnetic Cu 
spin configuration (GGA[fm]); (4) the spin-dependent 
GGA+U [S5J with a ferromagnetic Cu spin configura- 
tion (GGA+U[fm]). For the GGA+U[fm] calculations 
we considered the around mean field (AMF) version [SB] 
with values of U and Jh for the Cu ions of 6 eV and 1 eV 
respectively. The lattice constants, which are assumed to 
be well determined from experiments, were kept fixed for 
the structure relaxations while the optimization of the 
relative atomic positions was constrained by the symme- 
try of the Prima space group. 

For Cs 2 CuCl 4 , we also considered two optimization 
schemes with an antiferromagnetic Cu spin configuration: 
GGA[afm] and GGA+U[afm] ({7 = 6 eV, J H = 1 eV). 
The spin arrangement in this antiferromagnetic configu- 
ration is shown in Fig.[T](d), where the Cs 2 CuCl 4 unit cell 
was doubled in the b direction along the Cu chains. In or- 
der to produce such an arrangement, the symmetry of the 
supercell was lowered to the space group P21/c, with two 
inequivalent Cu atoms. Our choice of this particular an- 
tiferromagnetic configuration is due to its resemblance to 
the experimentally observed 120° ground state configura- 
tion |57| . The considered antiferromagnetic configuration 
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FIG. 1: (Color online) (a) The Cs 2 CuCl 4 and Cs 2 CuBr 4 unit 
cell where the CuA 4 tetrahedra (X = CI, Br) are highlighted, 
(b) Schematic triangular lattice of Cu spins. The gray lattice 
is below the black lattice, (c) Neighboring CuX 4 tetrahedra. 
Labels 1, 2 and 3 of the X atoms denote the three inequivalent 
X atoms: X(l), X(2) and X(3). As indicated, angles 712, 
713 correspond to the angles X(l)-Cu-A(2), X(l)-Cu-X(3), 
etc. (d) Cu spin configuration of a 1 x 2 x 1 supercell of 
Cs2CuCl 4 adopted for the structural optimization within the 
GGA+U[afm] scheme. 



3 



is collinear, which is beneficial in terms of computational 
effort, and fulfills the requirement that the strongest cou- 
plings J are satisfied and the second strongest couplings 
J' are partially satisfied. The optimization of the relative 
atomic positions with the two antiferromagnetic schemes 
was constrained by the symmetry of the P21/c space 
group. 

Additionally, for both systems we completed the LDA 
series of structural optimizations with LDA[fm] and 
LDA+U[fm] optimizations. However, due to the analo- 
gous behavior of the structural properties of CS2CUCI4 
and CS2CUBT4 observed within this series with those 
within the GGA series, the detailed analysis of the LDA 
series will be omitted here. 

The DFT structural optimizations were performed 
with the full-potential local-orbital (FPLO) code [rJHIBT] 
in the scalar relativistic approximation with up to 512 
k-points in the full Brillouin zone. 
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(b) Band gap 

Before presenting the structure optimizations with 
the different exchange-correlation functionals, we shall 
consider the experimental crystal structure and dis- 
cuss the electronic properties obtained with the dif- 
ferent exchange-correlation functionals. Calculations 
with the GGA[nm], GGA[fm] or GGA+U[fm] exchange- 
correlation potentials result in the electronic structure 
being either gapless (GGA[nm]) or gapped (GGA[fm] and 
GGA+U[fm]), as shown in Fig. [2} Allowing for spin po- 
larization opens a gap in both CS2CUCI4 and Cs2CuBr4 
[though in the latter compound the GGA[fm] gap is 
rather small (~ 0.03 eV)]. Upon introducing the onsite 
Coulomb repulsion within the GGA+U[fm], the gaps in 
both systems increase considerably [Fig.[2](c) and (f)]. 

Spin-dependent exchange-correlation functionals pro- 
vide a better description of CS2CUCI4 and Cs2CuBr4 elec- 
tronic properties as these functionals correctly reproduce 
the experimentally observed insulating ground state of 
the compounds. Moreover, as it will be shown in the 
next sections, it turns out that the presence of a band 
gap is important for the accurate determination of the 
equilibrium crystal structures. 



(c) Structural analysis 

We proceed now with the structural optimization. The 
crystal structures of CS2CUCI4 and Cs2CuBr 4 resulting 
after optimization with each of the schemes described 
above are presented in Appendix A. In our analysis of 
these structures we focus on geometry variations of the 
CUA4 (X — CI, Br) tetrahedron, which determine the 
strength of important exchange couplings. 



FIG. 2: (Color online) Total densities of states for major- 
ity and minority electron spin for (a) CS2CUCI4, calculated 
with GGA[nm], (b) Cs 2 CuCl 4 , calculated with GGA [fin], 
(c) Cs 2 CuCl 4 , calculated with GGA+U[fm], (d) Cs 2 CuBr 4 , 
calculated with GGA[nm], (e) Cs2CuBr4, calculated with 
GGA[fm], and (f) Cs 2 CuBr 4 , calculated with GGA+U[fm]. 
For an easier comparison, the two DOS's are plotted with op- 
posite signs. The calculations are performed with the LAPW 
code Wien2k [see Section II (d)]. In the case of GGA+U[fm], 
the around mean field scheme is employed, with U — 6 eV 
and J H = 1 eV. The Fermi level is set to zero. 



In CS2CUCI4 and Cs 2 CuBr4, the CuA 4 tetrahedron is 
distorted due to the Jahn- Teller effect and also due to 
the steric pressure originating from Cs + ions 0j. The 
Jahn- Teller effect results in a squeezing of the tetrahe- 
dron such that the X-Cu-X bond angles 712 and 733 in- 
crease and the X-Cu-X bond angles 713 and 723 decrease 
[Fig. [l](c)]. The steric pressure by Cs(2) on X(3) causes 
additional symmetry lowering by increasing 713 and de- 
creasing 733. In order to compare the strengths of the two 
types of distortions in different structures, one can define 
the Jahn- Teller deviation A JT as the difference between 
averages 7,(712 +733) and 5(713+723), 
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and the steric pressure deviations (5^ tcrlc and Sf' 

cstcric 2l7l2 — 733 I 
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As can be seen in Appendix A, the CS2CUCI4 struc- 
tures relaxed with GGA[fm] and GGA[afm] are very 
close as well as the structures relaxed with GGA+U[fm] 
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FIG. 3: (Color online) (a) The average Cu-X (X = CI, Br) bond distance dc u -x, (b) Jahn-Teller deviation A JT , (c) steric 
pressure deviation 5f crlc , (d) steric pressure deviation (5| tcrlc and (e) Cu-X-X angle a in the J superexchange bridge in the 
experimental as well as relaxed CS2CUCI4 and Cs2CuBr4 crystal structures. 



and GGA+U[afm]. The tetrahedron parameters in these 
structures are also similarly close. This demonstrates 
that, within the spin-dependent GGA and GGA+U, in- 
teratomic forces in CS2CUCI4 are very weakly dependent 
on the actual Cu spin configuration and in the follow- 
ing we will consider only ferromagnetic spin-dependent 
relaxation schemes for CS2CUCI4 and Cs2CuBr4. This 
result is a consequence of the rather small energy scale of 
magnetic interactions between Cu spins, which is much 
smaller than the difference between total energies usually 
involved in structural relaxations. 

In Tables [T| and pi} we present the tetrahedron parame- 
ters defined above for the experimental and relaxed crys- 
tal structures of CS2CUCI4 and Cs2CuBi'4, respectively. 
Some of them are additionally shown in Fig. [3] in or- 
der to facilitate the comparison of the different crystal 
structures. In both compounds, the crystal structures 
of the GGA relaxation series (GGA+U[fm], GGAffm] 
and GGA[nm]) are featured by a continuous variation 
of the X-Cu-X angles. In terms of these angles, the re- 
laxed crystal structures closest to the experimental struc- 
tures are the GGA+U [fm] structures for both CS2CUCI4 
and Cs2CuBi'4. However, in terms of Cu-X bond dis- 
tances, represented here as the averaged distance dc u -x, 
the GGA crystal structures are further away from the 
experimental structures than the LDA structures [see 
Fig. [3] (a)]. In the LDA structures (including the LDA[fm] 
and LDA+U[fm] ones, which are not listed in the tables), 
d-Cu-x are smaller than the corresponding parameters 
of the experimental structures, but the difference is less 
compared to the GGA series. Thus, the GGA relaxation 
tends to increase the bond distances in the CUX4 tetra- 
hedron while the LDA relaxation decreases them. 

In the GGA series, the Jahn-Teller-like distortion of 
the tetrahedron, which we evaluate through A JT , is 



TABLE I: Tetrahedron parameters for the CS2CUCI4 struc- 
tures and corresponding values of the Cu magnetic moment 
during structural relaxation. The angles are given in degrees, 
the deviations Sf oric 
Cu-Cl distance dcu 



and 82 erlc in percent and the averaged 
ci in An gstrom. 



exp GGA+U [fm] GGA[fm] GGA[nm] LDA[nm] 



712 
7i3 

723 
733 

A 



IT 



sf 



131.33 
101.67 

99.58 
126.79 

28.44 
153.99 

1.76 
1.04 



rfcu-ci 2.232 



130.97 
101.60 

99.43 
127.78 

28.86 
153.89 

1.23 
1.08 

2.263 

0.78^ s 



133.05 
100.67 

98.55 
130.41 

32.12 
155.21 

1.00 
1.06 

2.263 

0.50/xs 



133.66 
99.69 
97.86 

134.34 
35.22 

157.17 

0.25 
0.93 

2.267 





135.56 
98.26 
97.11 

138.56 
39.38 

159.28 

1.10 
0.59 

2.227 





strongest in the GGA[nm] crystal structures of both 
compounds and decreases monotonically as the Cu 
magnetic moment in the spin-resolved calculations in- 
creases (GGA[fm], GGA+U [fm]) [see Fig. [| (b)]. The 
Cu magnetic moment in CS2CUCI4 calculated with ei- 
ther GGA[fm] or GGA[afm] schemes is 0.50/1 s and in- 
creases to 0.78/is when the onsite Coulomb interaction is 
switched on within the GGA+U functional (Table [I]). In 
Cs2CuBr4, the corresponding values of the Cu magnetic 
moment are 0.42^ B and 0.73/xs [SS] (Table [n]). In the 
LDA[nm] structures, A JT takes the largest value. The 
distortion due to steric pressure, which is characterized 
by Sf CIlc and ^f " [Figs. [3] (c) and (d)], does not appear 
to follow any general rule. 

For the further discussion of the electronic structure it 
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TABLE II: Tetrahedron parameters for the Cs2CuBr4 struc- 
tures and corresponding values of the Cu magnetic moment 
during structural relaxation. The angles are given in degrees, 
the deviations <5f terlc and <5| torlc in percent and the averaged 
Cu-Br distance dcu-Br in An gstrom. 

~~ exp GGA+U[fm] GGA[fm] GGA[nm] LDA[nm] 



712 
713 
723 

733 
A JT 



6? 



130.40 
102.16 

99.93 
126.42 

27.52 
153.21 

1.55 
0.96 



dcu-Br 2.376 



130.06 
102.16 

99.75 
126.42 

27.28 
153.22 

1.42 
1.19 

2.407 

0.73/iB 



131.77 
101.70 

98.89 
128.05 

29.62 
154.03 

1.43 
1.40 

2.407 

0.42/is 



132.53 
100.48 

97.94 
133.07 

33.59 
156.53 

0.20 
1.28 

2.411 





133.72 
99.03 
97.34 

137.48 
37.41 

158.74 

1.38 
0.86 

2.368 





is useful to consider an additional structural parameter, 
namely, the angle in the J superexchange bridge Cu-X- 
X-Cu, shown in Fig. [11(c) as a. This angle is closely re- 
lated to 733 and behaves analogously in the various struc- 
tures [see Fig.[3](e) and Tab.|l|. The large variation of the 
superexchange coupling J in different CS2CUCI4 relaxed 
crystal structures is mainly attributed to the variation of 
a, as will be discussed in Section III. 
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(d) Electronic structure 

In this section, we present a detailed analysis of elec- 
tronic properties of CS2CUCI4 and Cs2CuBr4 for the 
experimental and relaxed crystal structures introduced 
in the previous section. We show the results obtained 
with the spin-independent GGA exchange-correlation 
functional as these are also used for the tight-binding 
parametrization. Calculations were performed with both 
the FPLO code as well as the linearized augmented plane 
wave (LAPW) scheme, as implemented in the Wien2k 
code [59]. The calculations with both codes are in good 
agreement and the data presented here are obtained with 
the LAPW code. 

Fig. [4] displays the density of states (DOS) for 
CS2CUCI4 and Cs2CuBr4 obtained using the experimen- 
tal structure. In both compounds, the hybridized Cu 3d 
and CI 3p/Bi Ap bands occupy the energy range between 
about -4 eV and eV [see, for instance, the atomic DOS 
in Figs. [4] (a) and (c)]. At the Fermi level, the Cu 3d xy 
states are half-filled as Cu is in a 3d 9 configuration. There 
are four Cu 3d xy bands in the bandstructure, which cor- 
responds to the number of Cu atoms per unit cell. A gap 
of approximately 4 eV separates the Cu and X (X = CI, 
Br) bands from the next unoccupied states [not shown in 
Figs. [4] (a) and (c)], which have significant Cs contribu- 
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FIG. 4: (Color online) Atomic species resolved densities of 
states (DOS) calculated from the experimental crystal struc- 
tures for (a) CS2CUCI4 and (c) Cs 2 CuBr 4 and the orbital pro- 
jected densities of states of Cu 3d for (b) CS2CUCI4 and (d) 
Cs2CuBr4. Energy is measured relative to the Fermi level _Ef. 



tion. 

Almost no contribution from Cs atoms to the DOS 
near the Fermi level is observed indicating a negligible 
hybridization of Cu with Cs. In particular, this indicates 
that exchange coupling J along the Cu chains in the b 
direction arises from the Cu-ATW-Cu hybridization. 

The Cu and X band manifold is an assembly of bond- 
ing states in the interval between -4 eV and -2 eV and 
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FIG. 5: (a) A schematic splitting of Cu 3d orbitals in the 
crystal field, (b) The path in the Cs 2 CuCl4/Cs 2 CuBr4 Bril- 
louin zone for bandstructure calculations (Figs. [6] and FFJ. (c) 
The local reference frame of a Cu atom, in which its orbital 
projected density of states is defined. 



antibonding states in the interval from -2 eV up to the 
Fermi level. The Cu antibonding states are split by the 
crystal field generated by X~ ions surrounding a Cu 2+ 
ion into the energetically lower Cu e a doublet (d x 2_ y 2 and 
d z z) and the energetically higher Cu ti g triplet (d xy , d xz 
and d yz ). Due to the Jahn- Teller-like uniaxial distortion 
of the tetrahedron, the t2 g triplet is further split into the 
degenerate d xz /d yz states and the half- filled d xy states. 
This splitting is schematically illustrated in Fig. [5] (a), 
and the orbital projected densities of Cu 3d states for 
the experimental CS2CUCI4 and Cs2CuBr 4 structures are 
presented in Figs. [4] (b) and (d), respectively. Note that 
the orbital designation is given according to the local 
reference frame of the CuA^ tetrahedron as shown in 
Fig. ©(c). 

In Figs. [6] and [7J we present the total DOS and 
bandstructures for the experimental as well as relaxed 
crystal structures (GGA+U[fm], GGAffm], GGA[nm], 
LDA[nm]) of CS2CUCI4 and Cs2CuBr4, respectively. In 
contrast to CS2CUCI4, the t^g states in Cs2CuBr4 are 
strongly hybridizing, which is indicated by the non- 
separable character of the overlap of the Cs2CuBr4 d xy 
and d xz /d yz bands in the bandstructure. 

We now focus on the comparison between the elec- 
tronic structures of CS2CUCI4 and Cs2CuBr4 calculated 
from different crystal structures. In both compounds, the 
choice of the relaxation scheme determines the degree 
of separation between the d xy and d xz /d yz bands and 
the dispersion of the d xy band. The separation between 
the d xy and d xz /d yz bands is reduced in the GGA[nm] 
crystal structure, compared to the LDA[nm] one, and 
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Energy (eV) 

FIG. 6: (Color online) DOS and bandstructures for various 
CS2CUCI4 structures (specified by the panel captions). In the 
bandstructure plots, the DFT calculated bands are in black 
(solid) lines and the tight-binding fits are in red (dashed) 
lines. The bandstructure path in the Brillouin zone is shown 
111 Fig.[5](b). Ener gy is again measured relative to the Fermi 
level E-p. 



keeps reducing within the GGA series as the value of the 
Cu magnetic moment, associated with a given relaxation 
functional, gets larger (GGA[fm], GGA+U[fm]). This 
trend is better seen in CS2CUCI4 where the d xy band is 
separated by a gap (except for the GGA+U[fm] struc- 



7 




-2.0 -1.5 -1.0 -0.5 0.0 Y r Z U X r T Y 
Energy (eV) 

FIG. 7: DOS and bandstructures of various Cs2CuBr4 struc- 
tures (specified by the panel captions). The bandstructure 
path in the Brillouin zone is shown in Fig. [5] (b). 



ture). We relate the variation of the degree of separation 
to the variation of the Jahn- Teller-like distortion of the 
CuX 4 tetrahedron. Our assumption is supported by the 
structural analysis in Tables [I] and |TT| as one observes there 
that the structures possessing more distorted tetrahedra 
demonstrate larger d xy and d xz /d yz band separation in 
the electronic structure. 

Analysis of the d xy band dispersion shows that the 



d xy bandwidth does not change significantly among the 
Cs2CuBr4 structures including the experimental one 
(Fig. [7|. In Cs 2 CuCi4, the d xy bandwidth of the exper- 
imental crystal structure is significantly narrower com- 
pared to that of the relaxed crystal structures, where it 
stays nearly uniform (Fig. [6| . The shape of the d xy bands 
varies in details from structure to structure in both com- 
pounds. This shape variation is accompanied by a vari- 
ation of the density of states distribution within the d xy 
orbital, which is easier to quantify. Thus, in both com- 
pounds, the weight of the d xy DOS gets shifted closer to 
the Fermi energy in the LDA[nm], GGA[nm], GGA[fm], 
GGA+U[fm] crystal structures sequence, which results 
in a more peaked appearance of the d xy DOS. The DOS 
at the Fermi level increases correspondingly. In the ex- 
perimental CS2CUCI4 structure, the DOS at the Fermi 
level is extremely sharp also due to the flatness of the 
d xy band. 



(e) TB model for Cs 2 CuCl 4 

The changes observed among the various electronic 
structures of CS2CUCI4 and Cs2CuBr 4 can also be dis- 
cussed on a quantitative level by mapping bandstructures 
to a tight-binding (TB) Hamiltonian: 



H 



TB 



E 



(3) 

where cj a and Cj i(T are, respectively, the creation and 
annihilation operators of electrons on site % and with spin 
<7, Uj denotes the hopping integrals and h. c. denotes 
Hcrmitian conjugate terms. 

Since in Cs 2 CuCi4 the Cu 3d xy bands at the Fermi level 
are well separated from the rest of the t 2ff manifold, its 
low-energy properties can be described by a one-band TB 
model. The hopping integrals are evaluated by Fourier- 
transforming the Hamiltonian Eq. ^ and mapping its 
eigenvalues to the four DFT Cu 3d xy bands. 

In Cs 2 CuBr 4 , the TB procedure is complicated by the 
hybridization of the overlapping d xy and d xz /d yz bands. 
In this three-band TB model is required. Since the 

three-band TB model includes new electronic degrees of 
freedom due to electron hopping between different types 
of orbitals, the number of independent hopping integrals 
in this model is considerably larger, compared to the one- 
band case. With an increasing number of model parame- 
ters, the TB parametrization by means of fitting becomes 
less reliable, and therefore we do not apply this method 
to Cs2CuBr4 in the present work. The rest of this sec- 
tion will be dedicated to the discussion of the CS2CUCI4 
single-band TB models. 

Though it was experimentally established [TT] that the 
spin interactions in Cs 2 CuCl 4 are predominantly within 



TABLE III: The TB model parameters in meV for the 
CS2CUCI4 Cu 3d xy band, calculated from the various 
CS2CUCI4 crystal structures. The hopping integral index cor- 
responds to the order of the neighbor. The less important 
interaction pathways, not present in Fig. [T](b), are shown in 
Fig. [9] (a) (Appendix B) . 
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the Cu layer along the J and J' paths [see Fig. [T] (b)], 
we observe that the electronic behavior modeled by 
the TB Hamiltonian show non-negligible interlayer hop- 
ping terms. We find that the d xy bands of all the 
Cs 2 CuCl4 structures can only be satisfactorily described 
with a minimal model that includes five hopping inte- 
grals. Three of them are the intralayer hoppings t and 
if and the interlayer hopping t" (t, t' , t" correspond to 
the interaction paths J, J', J", etc.), which have been 
considered in previous studies in the framework of spin 
Hamiltonians. The two new hopping parameters are the 
interlayer t\ and £3, also shown in Fig. [l] (b), their in- 
dices indicating the order of the interacting Cu neighbor. 
Note that t, t' and t" as well as J, J' and J" corre- 
spond to fifth, fourth and second nearest neighbour in 
the structure, respectively, while for all other ij and Ji 
the index corresponds to the order of the neighbour. In 
Table |III[ we present the values of relevant hopping in- 
tegrals for the various CS2CUCI4 structures. We observe 
that t\ and t% are in many cases comparable in magnitude 
to the hopping integral t' . We also note that the pres- 
ence of non-negligible interlayer couplings is revealed by a 
significant dispersion of the d xy bands along the a* direc- 
tion in fc-space (the X — T path in Fig. [6| . Additionally, 
we confirmed the importance of t\ and £3 by perform- 
ing NMTO downfolding calculations [SUj for some of the 
relaxed crystal structures. 

The total number of TB hopping integrals considered 
for an accurate description of the DFT bandstructure 
of CS2CUCI4 amounts to 11, plus the onsite energy /j, 
(Table III). Fig. [6] displays in dashed red lines the TB 
fit to the DFT Cu 3d xy bands, which are in solid black 
lines. Interaction pathways for the first seven hopping 
integrals are shown schematically in Figs.[T](b) and [9] (a) 
(in Appendix A). 

Comparing the TB models for different CS2CUCI4 



structures, we observe that the dominant hopping inte- 
gral t has the smallest absolute value for the experimental 
crystal structure and increases in the GGA crystal struc- 
ture series with a decrease of the Cu magnetic moment 
during relaxation (GGA+U[fm], GGA[fm], GGA[nm]). 
It has a maximum value in the LDA[nm] structure. At 
the same time, the second important coupling tf does not 
change significantly among the relaxed structures, but is 
by a factor of two smaller in the experimental crystal 
structure. The interlayer coupling i 3 is stable for all the 
structures, while t" and t\ vary considerably. 

We associate the variation of the hopping integral t in 
different CS2CUCI4 crystal structures with the variation 
of the CuCi4 tetrahedron geometry and, in particular, 
with the variation of the angle a in the Cu-Cl(3)-Cl(3)- 
Cu interaction path. In the relaxed CS2CUCI4 crystal 
structures, larger t values correspond to larger values 
of a and larger tetrahedron distortions. However, the 
experimental structure does not follow this rule; while 
the angles of the experimental and GGA+U[fm] relaxed 
structures are very similar, the t values differ by a factor 
of two. The reason why the experimental structure de- 
viates from the relaxed structures might be the influence 
of other tetrahedron parameters, such as details of the 
Cu-Cl bond distances. 

The observed relation between the dominant hopping 
integral t and the angle a is reasonable since a is the 
defining angle for the Cu 3d xy -Cl 3p-Cl 3p-Cu 3d xy hy- 
bridization. By considering perturbation theory on the 
onsite Coulomb repulsion U up to the second order, the 
effective Cu-Cu superexchange coupling can be obtained 
from t as J = Then, the relation between t and a 
fulfills the Kanamori-Goodenough rule stating 
that J reaches a maximum when the cation-anion-cation 
angle equals 180°. In the present case of the cation-anion- 
anion-cation (Cu-Cl-Cl-Cu) bridge, the four atoms get 
aligned along a straight line when a increases. 



III. EXCHANGE INTEGRALS 
(a) Computational details 

The exchange coupling integrals of the Heisenberg 
Hamiltonian, 



#H = 



E 

<i,j> 



Jij $i $ j j 



(4) 



can be obtained by means of DFT spin-resolved to- 
tal energy calculations [64 . Considering the differences 
between the energies of the ferromagnetic configura- 
tion and various antiferromagnetic spin configurations, 
B FM — E^ M , one derives a set of coupled equations for 
the couplings Jy. Following Eq. Q, antiferromagnetic 
exchange corresponds to positive Jij. 
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(b) Cs 2 CuCl 4 
Exp. structure 
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GGA+U[fm] structure 
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FIG. 8: (Color online) The spin exchange coupling constants as functions of 1/U for (a) the CS2CUCI4 experimental structure 
calculated with the atomic limit version of the GGA+U, (b) the CS2CUCI4 experimental structure (c) the CS2CUCI4 GGA+U[fm] 
relaxed structure, (d) the CS2CUCI4 GGA[fm] relaxed structure, (e) the CS2CUCI4 GGA[nm] relaxed structure, (f) the Cs2CuBi'4 
experimental structure, (g) the Cs2CuBr4 GGA+U[fm] relaxed structure, (h) the Cs2CuBr4 GGA[fm] relaxed structure, (i) the 
Cs2CuBr4 GGA[nm] relaxed structure. When not specified otherwise, the exchange couplings are obtained with the around 
mean field version of the GGA+U. The three sets of exchange couplings correspond to U = 8, 6 and 4 eV. Dashed lines mark 
the experimentally determined values of J, J' and J". 



For CS2CUCI4 and Cs2CuBr4, we take into account 
seven important couplings corresponding to the first 
seven hopping integrals in Table |III| which are J, J', 
J^' ff = J" + J 6 , Ji, J 3 , J 7 and J 8 [see Figs.[l](b) and[9](a)]. 
This choice defines the number of antiferromagnetic con- 
figurations to be calculated. The combined coupling J" s 
is introduced because the considered unit cell (which is a 
2x2x1 supercell) does not allow a separate calculation 
of the couplings J" and Jq but only their calculation as 
a sum. Since J" is presumably larger than Jq, J' e ' s gives 
an approximate value of J". 

The choice of the supercell is dictated by the peculiari- 
ties of the CS2CUCI4 and Cs2CuBr 4 magnetic sublattices. 
First, since the candidate for the largest coupling J con- 
nects Cu atoms that belong to adjacent primitive unit 
cells in CS2CUCI4 and Cs2CuBr4, the primitive unit cell 
has to be doubled in the b direction. Otherwise, J would 
be always canceled in any E FM — ijA™ difference. Also, 
in order to be able to discern the inequivalent couplings 
Ji and J3, we double the unit cell once more in the a di- 



rection and thus end up with a 2 x 2 x 1 supercell. In the 
supercell, we set eight out of 16 Cu atoms inequivalent 
in order to be able to arrange the required seven an- 
tiferromagnetic configurations within the same unit cell 
space group, which is P-l. It is important to stay within 
the same space group during total energy calculations 
for CS2CUCI4 and Cs2CuBr4 as in these compounds the 
exchange couplings are small and integration over differ- 
ently sampled Brillouin zones can affect the accuracy of 
the results. 

The seven antiferromagnetic spin configurations, to- 
gether with the set of coupled equations for the exchange 
couplings, are presented in Appendix C. The total en- 
ergy calculations were performed with the FPLO code. 
Test calculations with Wien2k confirm the results. In the 
FPLO code, we chose a 5 x 4 x 3 mesh of fc-points for 
the supercell Brillouin zone integration and kept other 
settings at default. The scheme to compute magnetic 
exchange for a given structure consisted of a series of 
total energy calculations within the AMF version of the 



10 



GGA+U, with U = 4, 6 and 8 eV and J H = 1 eV in 
all cases. This scheme was applied to the experimental 
structures of CS2CUCI4 and Cs 2 CuBr4 and also to the 
GGA[nm], GGA[fm] and GGA+U[fm] relaxed structures 
of the two compounds. 

(b) Results 

The calculated exchange couplings are plotted in 
Fig. [8j as a function of 1/U. For an easy comparison, 
we also mark by horizontal dashed lines of corresponding 
colors the experimentally determined values of J, J' and 
J" for Cs 2 CuCl 4 PU and of J and J' for Cs 2 CuBr 4 [TSj. 
Additionally, the exchange coupling values are provided 
in Appendix D. 

We first compare the Heisenberg models of CS2CUCI4, 
derived for the experimental as well as relaxed crystal 
structures [Figs. [|] (a)-(e)]. The common feature for 
all these structures is that their spin models are two- 
dimensional, with J being the leading interaction fol- 
lowed by J' . 

For the experimental structure [Fig. [8] (b)], the ra- 
tio J' / J is similar to the experimentally determined ra- 
tio [HI [65], while the absolute values of J and J', which 
are calculated in the present work, are off by a factor of 
almost 3. In this structure, the two-dimensionality is less 
pronounced due to the interlayer antiferromagnetic inter- 
action J3 and the ferromagnetic interaction J7, which are 
competing with J'. The GGA[nm] structure [Fig. [8] (e)], 
on the other hand, tends to behave as a ID rather than 
2D system as here the J'/J ratio is 0.15 at U = 6 eV. We 
assume that U = 6 eV should be close to the U value for 
the Cu 3d electrons in the LAPW GGA+U scheme (see, 
for instance, Ref . l6"6"lwhere the authors successfully repro- 
duce the electric field gradients at Cu 2+ ions in a num- 
ber of Cu 2+ oxides and halides by performing LDA+U 
calculations with U cS = U - J = 5 eV). The J'/J ra- 
tios at U = 6 eV in the experimental, GGA+U [fm] and 
GGA[fm] structures are, respectively, 0.30, 0.47 and 0.27. 
These ratios vary slightly for other U values. 

The differences in the J'/J ratios among the relaxed 
structures are mainly due to a strong variation of the 
coupling J, which adopts the values 0.384, 0.568, 0.932 
meV for U = 6 eV in the GGA+U [fm], GGA[fm] and 
GGA[nm] structures, respectively, whereas J' decreases 
only slightly for the same sequence of structures. 

Judging by the proximity of the theoretically derived 
Cs 2 CuCl4 Heisenberg model (this work) to the experi- 
mentally determined one [llj . the calculations with the 
GGA[fm] and GGA+U[fm] optimized structures provide 
the most satisfactory results. In both cases, the abso- 
lute values of J and J' as well as their ratios are close to 
experiment for typical values of U between 6 and 8 eV 
(for GGA+U[fm] U « 6 eV yields the best agreement, 
while in the case of GGA[fm] interpolation of the results 



between 6 and 8 eV results in almost perfect agreement 
with experiment for U ~ 7 eV). Moreover we find non- 
negligible interlayer couplings as shown in Appendix D 
(bold- face values). 

For the spin models calculated with the relaxed 
CS2CUCI4 structures, the ratios of antiferromagnetic ex- 
change couplings are in good agreement with the ratios 
of the corresponding squared hopping integrals. This is 
an indication that ferromagnetic contributions to these 
couplings are small. 

In contrast to CS2CUCI4, the Cs2CuBr4 effective spin 
model derived using the experimental crystal structure 
[Fig. [8] (f )] agrees well with the experimentally estimated 
model parameters in the interval of U values around 
6 eV. We obtain that Cs2CuBr4 is a two-dimensional 
system, with dominant antiferromagnetic couplings J 
and J' and considerably smaller interlayer couplings. At 
[7 = 6 eV, the obtained J'/J ratio equals 0.64, which 
compares well with the experimental result of 0.74. Also, 
the model based on the GGA+U [fm] relaxed structure 
[Fig. [8] (g)] gives similar results for the exchange parame- 
ters. The similarity between these two models is another 
feature that distinguishes Cs2CuBr4 from CS2CUCI4. In 
the latter case, as seen above, the exchange couplings 
of the experimental model [Fig. [8] (b)] are considerably 
smaller than those of the GGA+U [fm] relaxation model 
[Fig- 1(c)]. 

Overall, the Heisenberg models for Cs 2 CuBr4 obtained 
within the GGA relaxation series [Figs. [8] (g)-(i)] follow 
the same behavior as the corresponding CS2CUCI4 mod- 
els, which is characterized by increasing J and decreasing 
J'/J in structures that have been relaxed with a smaller 
Cu magnetic moment. 

To conclude this section, we briefly comment on the 
performance of the atomic limit (AL) version [54] of the 
GGA+U exchange-correlation functional, which is an al- 
ternative to the around mean field version. The exchange 
couplings of the experimental Cs 2 CuCl4 structure cal- 
culated with the AL double counting correction differ 
considerably from those obtained with the AMF dou- 
ble counting correction [compare panels (b) and (c) of 
Fig. [8] . The AL calculated exchange couplings J and J' 
do not behave linearly with 1/U, as expected from the 
J = At 2 /U relation, valid for these antiferromagnetic cou- 
plings. Therefore, we restrict ourselves in the remaining 
discussion to calculations with the AMF double counting 
correction. 



IV. DISCUSSION 

From our analysis of the structural and electronic prop- 
erties of CS2CUCI4 and Cs2CuBr4 we obtained the follow- 
ing results: Within the sequence LDA[nm], GGA[nm], 
GGA[fm] and GGA+U[fm] of functionals used for the 
structural optimization, the Jahn-Teller-like distortion of 
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the CuXt tetrahedron is reduced, which is accompanied 
by an increase of the Cu-X(3)-A(3) angle a in the J 
superexchange path. These structural changes lead to 
a considerable variation of the exchange coupling J and 
the corresponding hopping integral t. 

(a) Importance of magnetism and electronic 
correlations 

We find that the choice of the functional used in the 
structure optimization is crucial for the correct modeling 
of the properties of these two compounds, especially in 
the case of CS2CUCI4. First of all, opening of a gap at 
the Fermi level by introducing various magnetic struc- 
tures seems to be a necessary ingredient in the structure 
optimization, as can be seen from the relatively bad per- 
formance of the LDA[nm] and GGA[nm] optimizations. 
Since in these schemes the system is gapless, which is in 
stark contrast to the true ground state of CS2CUCI4 ob- 
served experimentally, the strong instability at the Fermi 
level (visible in the peak in the DOS) has to be par- 
tially relieved by forcing an unphysically strong struc- 
tural distortion, leading to an improper determination 
of lattice structure. Furthermore, a change of the size 
of the gap through onsite correlations and corresponding 
localization of the Cu magnetic moments also influences 
the results of the structure optimization. Although one 
might have expected a better performance of the spin- 
dependent optimization schemes (due to the magnetic 
nature of the compounds), the strong impact on the cal- 
culated exchange couplings is rather unexpected. The 
reason is presumably the small values of the exchange 
constants. 



(b) Comments on the quality of the experimentally 
determined CS2CUCI4 and Cs2CuBr4 crystal 
structures 

Comparison of electronic behavior and effective mod- 
els of the various relaxed Cs 2 CuCi4 structures with those 
of the experimental structure indicates that the experi- 
mental determination of the CS2CUCI4 crystal structure 
[3] was probably not sufficiently accurate, at least, as 
far as calculations of microscopic models are concerned. 
This would explain the huge differences between our de- 
rived CS2CUCI4 Heisenberg model when the experimental 
crystal structure is considered, and the model Coldea et 
al. [11] obtained from fitting to neutron data and which 
was corroborated by a number of studies [lOl HH [27J [69] . 

Accurate determination of the CS2CUCI4 crystal struc- 
ture might be complicated due to the presence of non- 
stoichiometric hydrogen containing compounds, presum- 
ably HC1 or H0 2 , detected in this material and not de- 
tected in Cs 2 CuBr 4 [SHE]- 



The experimentally determined Cs2CuBr4 crystal 
structure, on the other hand, is accurate enough. We 
suggest therefore that Cs2CuBr4 can be regarded as a 
reference system for choosing the relaxation scheme, suit- 
able also for describing CS2CUCI4. We conclude from the 
data analysis and from the physical considerations that 
such a relaxation scheme is the GGA[fm](GGA+U[fm]), 
with < U < 6 eV. 

(c) Comparison with the experimentally determined 
microscopic models 

In view of the arguments presented above, one should 
refer to the results obtained with either GGA[fm] or 
GGA+U[fm] relaxed structures when discussing the real- 
istic spin model for CS2CUCI4. Our calculations confirm 
for both systems that their spin models are 2D, with 
the intraplane J and J' being the leading interactions, 
and Cs2CuBr4 showing a higher degree of frustration 
than CS2CUCI4. The model from Coldea et al. [11] . with 
J = 0.374 eV and J'/ J = 0.34, is close to our DFT 
models derived with the GGA[fm] and GGA+U[fm] re- 
laxed structures and remains a valid model for CS2CUCI4. 
In Cs2CuBr4, where the experimental crystal structure 
seems to be much more reliable, we can propose an 
approximate DFT model, obtained as a generalization 
of the models of Figs. [8] (f) and (g) for U « 6 eV: 
0-8 < J < 0.9, 0.5 < J' I J < 0.65. This model is quite 
close to the model by Ono et al. [15], with J = 0.97 eV 
and J' I J = 0.74. 

Additionally, in view of the recent theoretical studies 
by Starykh et al. [37] IHH] > who demonstrated the impor- 
tant role of the relatively weak Dzyaloshinskii-Moriya in- 
teraction and interlayer coupling in determining the mag- 
netic behavior of CS2CUCI4, our results for the effective 
Cs 2 CuCl4 model reveal a possible relevance of a number 
of interlayer couplings, besides J", which are of compa- 
rable strength (Ji, J3 and J7). Also, we find the next- 
nearest chain neighbor hopping integral tu to be quite 
large, thus supporting another suggestion by these au- 
thors. 

Finally, it is important to emphasize that the electronic 
TB models for the two compounds involve sizable inter- 
layer hoppings, comparable with the intralayer ones, so 
that in terms of electronic degrees of freedom the two 
compounds have 3D behavior. In the case of Cs2CuBr4, 
this can be seen from the large dispersion of the bands 
along the a* direction in fc-space, even though we didn't 
derive a TB model for this material. 

(d) Structural relaxation 

Finally, we comment on our decision to perform struc- 
tural relaxations with fixed lattice constants and using 
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the room temperature data. It is generally known that 
the GGA tends to overestimate the unit cell volume while 
the LDA underestimates it. Therefore, we preferred to 
rely on the experimental lattice constants, which are usu- 
ally determined with high accuracy. We performed a test 
GGA[nm] relaxation of the CS2CUCI4 crystal structure 
with lattice constants measured at 0.3 K [TT] and found 
that the electronic structure and the TB model of the 
resulting CS2CUCI4 structure are close to the electronic 
structure and the TB model of the GGA[nm] Cs 2 CuCl4 
structure, relaxed with the room temperature lattice con- 
stants. Therefore, all calculations have been performed 
with the room temperature lattice constants for both 
compounds. 

V. SUMMARY AND CONCLUSIONS 

In summary, we have performed DFT calculations for 
CS2CUCI4 and Cs2CuBr 4 . Our study shows that the ex- 
change coupling constants of these compounds exhibit a 
strong dependence on subtle details of the crystal struc- 
ture, especially on the geometry of the CuXi tetrahcdra. 
Depending on the structural model, we observe a large 
variation of the derived exchange couplings and their ra- 
tios, resulting in completely different spin model Hamil- 
tonians. One reason for this unusual sensitivity are the 
fairly small absolute values of the exchange couplings, 
with the largest coupling constant J being below 5 K. 

One important motivation for our detailed study is the 
failure of the experimental structure published in Rcf. 3 
to correctly describe the magnetic behavior of Cs 2 CuCl 4 . 
Calculations with the experimental structure provide too 
small exchange coupling constants with fairly strong in- 
terlayer couplings in contrast to the pronounced 2D char- 
acter observed in experiment. This indicates (together 
with the fairly large forces acting on the atomic positions) 
that a better characterization of the CS2CUCI4 crystal 
structure is necessary. Only after structural optimiza- 
tion with spin-dependent GGA and GGA+U functionals, 
we obtain an overall good agreement with the exchange 
couplings obtained from experiment [llj . 

In contrast, the leading exchange couplings for 
Cs2CuBr4 obtained from our calculations are in good 
agreement with those derived from experiment, indepen- 
dently of whether we use the experimental structure or 
structures from spin-resolved optimizations. 
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APPENDIX A: STRUCTURAL PARAMETERS 
OF CS2CUCL4 AND CS2CUBR4 OBTAINED 
WITHIN DIFFERENT OPTIMIZATION 
SCHEMES 

Below, we provide the Cs 2 CuCi4 and Cs 2 CuBr 4 rela- 
tive atomic positions obtained after structural optimiza- 
tion within different schemes. For a quick reference, we 
also cite Ref. [3] for the experimentally found structure of 
Cs 2 CuCl 4 and Ref. H for that of Cs 2 CuBr 4 . 

In the case of CS2CUCI4, for the experimental struc- 
ture and structures relaxed with non-spin-resolved and 
ferromagnetic calculations, the lattice constants are a = 
9.769 A, b = 7.607 A, c = 12.381 A and the space group 
is Prima. The structures relaxed with antiferromagnetic 
calculations (GGA [arm] and GGA+U[afm]), for which 
the relaxation was constrained by the symmetry of the 
P21/c space group in a supercell, were found to even- 
tually belong to the same space group P21/c but in a 
reduced cell, with the same unit cell parameters as those 
of the original full-symmetry unit cell of the compound. 
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In the case of Cs2CuBr4, the lattice constants are a = 
10.195 A, b = 7.965 A, c = 12.936 A and the space group 
is Pnma. 
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APPENDIX B: FURTHER INTERACTION 
PATHWAYS IN CS 2 CUCL 4 CONSIDERED IN 
THE TB MODEL 



Ji 
Je 



The interaction pathways for spin exchange coupling 
constants J 7 , Je, Ju, Jig and J22 are shown in Fig. [9] (a). 



APPENDIX C: ANTIFERROMAGNETIC SPIN 
CONFIGURATIONS AND CORRESPONDING 
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Given that the Cu atoms in the 2x2x1 super- 
cell of CS2CUCI4 or Cs2CuBr4 are labeled as shown in 
Fig.[9](b), the seven antiferromagnetic (or ferrimagnetic) 
spin configurations, considered in order to calculate ex- 
change couplings J, J', J' c ' s , J\, J3, J 7 and J 8 by the 
total energy difference method, are the following: 
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Each of these spin configurations leads to a corresponding 
equation: 



FIG. 9: (a) Interaction pathways J7, Je, Ju, Jia and J22. 
(b) Labeling of the Cu atoms in the 2x2x1 supercell of 
CS2CUCI4 or Cs2CuBr 4 . The black and gray circles denote 
Cu atoms that belong to adjacent Cu layers, parallel to the 
be plane. Since the supercell contains two unit cells along the 
a axis, the black and gray Cu lattices are doubled, which is 
not visible in the be projection. The two encircled figures, one 
on top of the other, label the Cu atoms that have common 
y and z coordinates, but whose x coordinates differ by the 
lattice constant a such that the top figure refers to the Cu 
atom with larger x coordinate. 
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where E t = E kM - Ef bM , i = !,...,!, with E being 
the energy of the supercell in the ferromagnetic config- 
uration of Cu spins and E^ FM being the energy of the 
supercell in the antiferromagnetic configuration i. 



APPENDIX D: EXCHANGE COUPLINGS 



The exchange couplings provided below for the ex- 
perimental as well as the GGA[nm], GGA[fm] and 
GGA+U[fm] relaxed CS2CUCI4 and Cs2CuBr 4 crystal 
structures are given in meV. We have marked in bold face 
the parameter values that should provide a realistic de- 
scription of the spin models for CS2CUCI4 and Cs2CuBr4. 
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Cs 2 CuCl 4 : experimental structure [Fig. [5] (b)] 
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Cs 2 CuCl 4 : GGA+U[fm] relaxed structure [Fig.[8](c)] 
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Cs 2 CuCl 4 : GGA[fm] relaxed structure [Fig. [5] (d)] 
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Cs 2 CuCl 4 : GGA[nm] relaxed structure [Fig. [5] (e)] 
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Cs 2 CuBr 4 : experimental structure [Fig. [5] (f)] 
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Cs 2 CuBr 4 : GGA+U[fm] relaxed structure [Fig. [S] 
(g)] 
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Cs 2 CuBr 4 : GGA[fm] relaxed structure [Fig. [5] (h)] 

7=4 eV 7=6 eV 7=8 eV 



./ 


2.7229 


1.2325 


0.4076 


./' 


0.9234 


0.4382 


0.1682 


7" 


0.0554 


0.0219 


0.0122 


Jl 


0.1467 


0.0492 


0.0166 


Ja 


0.1710 


0.1101 


0.0468 


Js 


-0.0158 


-0.0005 


0.0004 


Ji 


-0.0394 


-0.0183 


-0.0089 



Cs 2 CuBr 4 : GGA[nm] relaxed structure [Fig. [|] (i)] 
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